Accelerated cosmic expansion in a scalar-field universe 
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We consider here a spherically symmetric but inhomogeneous universe filled with a massless 
scalar field. The model obeys two constraints. The first one is that the gradient of the scalar field 
is timelike everywhere. The second constraint is that the radial coordinate basis vector is a unit 
vector field in the comoving coordinate system. We find that the resultant dynamical solutions 
compose a one-parameter family of self-similar models which is known as the Roberts solution. The 
solutions are divided into three classes. The first class consists of solutions with only one spacelike 
singularity in the synchronous-comoving chart. The second class consists of solutions with two 
singularities which are null and spacelike, respectively. The third class consists of solutions with 
two spacelike singularities which correspond to the big bang and big crunch, respectively. We see 
that, in the first case, a comoving volume exponentially expands as in an inflationary period; the 
fluid elements are accelerated outwards form the symmetry center, even though the strong energy 
condition is satisfied. This behavior is very different from that observed in the homogeneous and 
isotropic universe in which the fluid elements would move outwards with deceleration, if the strong 
energy conditions are satisfied. We are thus able to achieve the accelerated expansion of the universe 
for the models considered here, without a need to violate the energy conditions. The cosmological 
features of the models are examined in some detail. 

PACS numbers: 04.20.Dw,04.20. Jb,04.70 Bw 
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I. INTRODUCTION 



The big bang universe model has succeeded to explain 
three prominent observational facts, i.e., the Hubble's 
law, the abundance of light elements, and the 2.7K cos- 
mic microwave background radiation (CMBR) lj. Fur- 
ther, within this model, the formations of various struc- 
tures in the universe, i.e., stars, galaxies, clusters of 
galaxies, etc, seem to be naturally explained if we as- 
sume the existence of dark matters In the original 
big bang model, there are three basic unnatural aspects, 
which arc known as the horizon problem, flatness prob- 
lem, and the monopolc problem. Now it is widely be- 
lieved that these three problems are the evidences for 
existence of the so called inflationary period in the past 
of our universe 0, 0, IH- In this period, by virtue of 
the accelerated cosmic volume expansion, initial inho- 
mogencities are stretched, the space is made sufficiently 
flat, and at the same time, moderate quantum instabil- 
ities generate the seeds of the present structures in the 
universe. The theoretical prediction of the inflationary 
scenario basically agrees well with the WMAP data on 
CMBR@. 

Recent observational results also imply that our uni- 
verse has again entered into another inflationary period 
at about z ~ 0.3. This is known as the dark energy 
problem. The WMAP data and the observations of the 
distance-redshift relations of supernovae (tI \§i. r9l. [Toj mean 
that the speed of the cosmic volume expansion is nec- 
essarily accelerated if the universe is homogeneous and 
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isotropic. Further, if general relativity is the correct the- 
ory of gravity, then the accelerated cosmic expansion im- 
plies the existence of a dark energy component which 
does not satisfy the strong energy condition, i.e., effec- 
tively it has to have an equation of state given as P = wp 
with w < —1/3. It is far from clear today what such 
a dark energy component could be. The simplest an- 
swer to this conundrum could be the existence of a non- 
vanishing cosmological constant in the universe, since in 
that case, only one constant can explain all of the obser- 
vational data. However, it does not at all appear to be 
easy to explain the origin of such a cosmological term. 
Therefore, the issue is under an active investigation and 
several approaches to resolve the problem are being pur- 
sued simultaneously. These include the existence of a 
quintessence field, using an alternative theory of grav- 
ity (sec, for example Ref . fill] ) . and also inhomogeneous 
universe models, as we shall discuss below. 

In this paper, we consider a self-similar model for the 
spherically symmetric but inhomogeneous universe, filled 
with a massless scalar field which is minimally coupled 
to gravity. This solution to Einstein equations has been 
discussed in the literature by Roberts fl2l . Brady [l3[ 
and Oshiro, Nakamura and Tomimatsu ,14|. However, 
these authors used this solution only to study the grav- 
itational collapse of a massless scalar field and used it 
as a model for the collapsing matter. In this paper we 
will focus on the cosmological significance of this solu- 
tion. All of the physically reasonable energy conditions 
hold for this model universe. Although this solution does 
not directly relate to our present universe, there are cer- 
tain interesting properties that the models exhibit as we 
shall point out here. In the region where the gradient of 
the scalar field is timelike, the system is equivalent to a 
osaM-rtt.pejjiiect fluid, given by the equation of state P = p, 
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and we can consider motions of the fluid elements. As we 
show here, in this solution the fluid elements are acceler- 
ated outwards, even though the strong energy condition 
holds. 

Here we should mention the models for a dust-filled 
inhomogeneous universe, which have been considered to 
resolve the dark energy problem without introducing any 
exotic matter-components, or without invoking any alter- 
native theories of gravity. Mainly two approaches have 
been used in this connection. The first one invokes the 
backreaction effect of inhomogeneities in the universe, 
which causes the acceleration of the cosmic volume ex- 
pansion in an average sense [IH EH, E3, EH, EH, Hfj, HH ■ In 
the other approach, one invokes a spherically symmetric 
but inhomogeneous velocity field which causes the ap- 
parent acceleration of the cosmic volume expansion [2 2L 
IM EI HI H, HE Hi- At present, it seems to be dif- 
ficult to resolve the dark energy problem using the first 
approach (see e.g. [2^, [3(| HH, H2] for a discussion on this 
issue), and also there are some arguments which discuss 
the viability or otherwise of the second approach [H, [H| . 

We would like to emphasize that the present solution 
that we discuss and consider here is different from both 
of the above possible alternatives that have been used 
for dust-filled inhomogeneous universe models. In the 
present solution, each fluid element is accelerated out- 
wards in the real and exact sense from the observer at 
the center of symmetry. 

The paper is organized as follows. In Sec. II, we briefly 
review the massless scalar field minimally coupled to 
gravity, and then, in Sec. Ill, the basic equations in a 
comoving reference frame are presented. In Sec. IV, by 
imposing a constraint that the radial coordinate basis 
vector is a unit vector field, we obtain the solution. We 
believe our method gives some insights into the proper- 
ties of the solution we derived this way. In Sec.V, we 
show that the charts of the solutions obtained in some 
of the classes we considered do not cover whole of the 
spacctime. Therefore, we perform two kinds of exten- 
sions respectively, in Sec. VI and Sec. VII. In Sec. VI, wc 
also discuss and analyze the global structure of the so- 
lutions. Finally, Sec. VIII is devoted to a summary and 
discussion. 

In this paper, we adopt the abstract index notation(37j; 
the Latin indices denote the type of the tensor, whereas 
the Greek indices mean the components with respect to 
the coordinate basis. We use the units 8ttG = c = 1. 



II. MASSLESS SCALAR FIELD 

The massless scalar field 4>(x a ) minimally coupled to 
gravity on a spacetime (M, g a b) is described by the La- 
grangian, 



The corresponding Eulcr-Lagrange equation is 4>;abg ah = 
0, and the stress-energy tensor for the scalar field, as 
calculated from this Lagrangian is, 



Tab = <t>-a<t>-b - 2 gab i ( l ) ^ ( t ) - d 9 cd ) 



(2) 



In case that the gradient of the massless scalar field c/>- a 
is timclikc or spacclike, we can introduce a unit vector 
field 



N 



(3) 



where N = ^J\4>' b 4>-b\- Then the stress-energy tensor of 
the massless scalar field is rewritten in the form 



N 

Tab = — {2u a u b + eg a b) ■ 



(4) 



where e = — 1 for spacelike (f>- a , while e = +1 for time- 
like <j). a - Clearly, u a is one of eigenvectors of T ab . In 
this case, the massless scalar field is categorized into a 
Type 7 matter field (38j. i.e., it admits one timelike and 
three spacelike eigenvectors. So, at each point q G M, 
we can then express the tensor T ab in terms of an or- 
thonormal basis (Eft^E®, E 2 ,Eft), where Eft is a timelike 
eigenvector with an eigenvalue p and E a A (A = 1,2,3) are 
three spacelike eigenvectors with eigenvalues pa- Here p 
represents the energy density of the scalar field as mea- 
sured by an observer with a 4-velocity Eft at q, and the 
eigenvalues pa represent the principal pressures in three 
spacelike directions E a A . 

In the case of spacelike </> ;a , we can put Ef — u a . Then 
the stress-energy tensor Q is written in the form 



T 



ab 



N 2 



{Eft Eft + EiEl - E a 2 E b 2 - EftEl) . (5) 



It is seen from the above form that p = p\ = N 2 /2 > 0, 
whereas p 2 = pz = —N 2 /2 < 0. 

In the case of timelike 4>- a , we can put Eft = u a . Then 
the stress-energy tensor @ is written in the form 

N 2 

T ah = — (EftEft + E a x E\ + E a 2 E\ + EftE b 3 ) . (6) 

It is seen from the above form that p = p\ = p 2 = p% = 
N 2 /2 > 0. This is equivalent to the stress-energy tensor 
of the stiff perfect fluid. 

In the case that 4>- a is null, the stress-energy tensor 
takes the following form, 



Tab = (rk a kb, 



(7) 



where k a is a null vector field and a is a non-negative 
function. This is equivalent to the stress-energy tensor 
of the null dust, and is categorized into a Type //matter 
field. 



III. BASIC EQUATIONS 



£ = - \$-A;bg 



ab 



(i) 



Hereafter, we focus on the spherically symmetric sys- 
tem. 
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A. Comoving coordinate system 

If the gradient of the scalar field is not null, as men- 
tioned in the previous section, the stress-energy tensor 
has four eigenvectors (Eq,E%). In this case, we can 
choose the spherically symmetric coordinates (t, r, 6, ip) 
such that the coordinate basis vectors are parallel to 
these eigenvectors. This is the comoving coordinate sys- 
tem in the sense that there is no energy flux, i.e., the 
time-space components of the stress energy-tensor are 
vanishing. The line element in this coordinate system is 
then written in the following form, 

ds 2 = _ e Mt,r) dt 2 + e m,r) dr 2 + R 2^ (g) 

where dfl 2 is the line element on a unit two-sphere. It 
is noted that there are still two scaling freedoms of one 
variable left in t and r. In the models considered here, R 
is a monotonically increasing function of r. 

In the spherically symmetric system, a non-trivial 
time-space component is 



(9) 



where (') denotes the partial derivative with respect to t 
and (') with respect to r. Thus, in the comoving system, 
one of the <f> and <f>' must vanish. 

Hereafter, as mentioned in Sec. I, we assume that </> ;a is 
everywhere timelike, i.e., <f> = 4>{t), and thus this system 
is equivalent to that of the stiff perfect fluid. Here it 
is worth noting that if the norm of 4>- a does not have 
a definite sign, the comoving coordinate system cannot 
cover whole of the spacetime. 

The components of the stress-energy tensor are then 
given by, 



T^u = diag[-p, p, p, p] , 



where 



p = — e 
H 2 



-2vi2 



(10) 

(11) 



F = --e- 2v <j) 2 R 2 R , 



(R 2 e^- V 4) = , 



2R' + R'^+R^- = 



(16) 



(17) 



(18) 



Here the function F = F(t, r) has an interpretation of 
the mass function for the collapsing cloud, and it gives 
the total mass in a shell of comoving radius r on any 
spacelike slice t = const. 



IV. SOLVING THE EINSTEIN EQUATIONS 

The function R(t, r) is the area radius of a shell la- 
beled V at an epoch 'f. For the sake of definiteness, 
let us consider the situation of an expanding universe, 
so we have R > because we are considering here the 
expanding branch of the solutions. If R changes sign 
then that corresponds to the recollapse of the field dur- 
ing the dynamical evolution of the universe. We note 
that equation (TTT|) is the Klein-Gordon equation for the 
scalar field, which is a part here of the Einstein equations 
via the Bianchi identities. We can integrate this equation 
to get 



= rV(r), 



(19) 



where f(r) is an arbitrary function of integration. So- 
lution of these equations, subject to the initial data and 
energy conditions, would determine the time evolution of 
the system. 

We now construct a class of solutions with ip = 0. This 
is the second assumption that we make here, namely that 
the radial coordinate basis vector is an unit vector field. 
From ([TBI), we have 



B. Einstein equations 

The dynamic evolution of the initial data, as specified 
on a spacclike surface of constant time is determined by 
the Einstein equations. For the metric (J5J) , using the 
definitions 



G(t,r) 
H(t,r) 



-2,J 



-2// 



(Rf 



>\2 



and 



F = R(l-G + H), 



(12) 
(13) 

(14) 



the independent Einstein equations for the massless 
scalar field are, 



F' 



1 



-2vl2 r>2 



R R' , 



(15) 



H = H(t) 



(20) 



Also, using Eqs. (fT9|) and (|20|) , Eq. (flB) can be integrated 
to give 



l r 4 / 2 (r) 



(21) 



where h(r) is an arbitrary function of integration. To- 
gether with the above equation, Eq. (fT5)) leads 



1_V , (r 2 /)' 1 , h' 

R. 



0. 



r 2 f R r 4 f 2 
The solution of the above equation is given by 

r 2 f(r) 



R 



g(t)+ P (r)' 



(22) 



(23) 
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where g(t) is an arbitrary function of integration, and 

h'(r) 



p(r) = - 



2 f(r) 



dr. 



Here, by a scaling freedom of t, we set 
0(t) - ipt, 



(24) 



(25) 



where ip is a constant with a dimension of length-inverse 
(note that the scalar field is dimcnsionlcss in our unit). 
Then we have (pe~ v R 2 = r 2 f(r) from Eq. (fT9|) . Further, 
using Eqs. (fll?)) and (J2DJ), we have 



^ff(i) 



j r 4 f 2 (r)' 
Using Eq. (f2"3")l , the above equation implies 
a 2 (i)=^ 2 ff(i). 



(26) 



(27) 



Substituting the above results into Eq. (fT4|) . we have 



I, , 2 ft 



2 r z / 



(r 2 /)'( 3 +p) + ft' 



{g + pf 



(28) 

The left hand side of the above equation depends on only 
t, and thus the right hand side should also depend on only 
t. This equation gives a constraint on the functions / and 
ft. 

As long as we focus on the dynamical situations, we 
can see from Eq. (|2"5|) that the most general allowed form 
of r 2 /(r) is 



rf(r) = a, 



(29) 



where a is a positive dimcnsionlcss constant, and further 
ft = (see Appendix Ell. Then, Eq. ([28]) leads to 



<p- x g{t) = ± 



1 9/ . a 2 

2 sit) + 7W, 



(30) 



We take the negative sign in the equation so that R > 0. 
This can be easily solved, and we have, 



(31) 



where C is an integration constant. So the metric we get 
is 




ds 1 



Also, we have 



{(par) 2 2 2 



F = —arg(t) and 



(or) 2 

1 g 4 (t) 

2 (ar) 2 



dn 2 



(32) 



(33) 



From Eq. (|3"3"j) . we see that the energy density p diverges 
at r = if g does not vanishes. Moreover, even for r > 0, 



p diverges if g diverges. Since the Ricci scalar is propor- 
tional to p , g = oo corresponds to the scalar polynomial 
singularity [38j. 

Here, we choose C so that C = \j2cP- — 1 for a 2 > 
1/2, C = 1 for a 2 = 1/2, and C = \[\ - 2a 2 for a 2 < 
1/2. Then, we have following three distinct classes of 
solutions. 



A. The case of a 2 > 1/2 



In this case, we obtain 



(t)= 1 - \]2a 2 - 1 sinh ( V2ipt 



The domain of time t is — oo < t < t\>, where 



*b 



2\/2> 



In 



y/2a+ 1 



V^a- 1 



(34) 



(35) 



For t — > — oo, g diverges, and thus the energy density and 
the Ricci scalar also diverge for r > 0. As will be shown 
in the next section, t — > — oo is not infinity, and hence this 
should be regarded as an initial singularity (big bang). 
At t — ib, g vanishes, and thus any 2-dimensional spheres 
with positive comoving radii have infinite area, since the 
area A with a comoving radius r is given by 



A = Air(ar) 2 / g 2 



(36) 



As will be shown later, t — is on the future null infinity 
except at r — 0; a "point" (t, r) = (ib, 0) is a sphere with 
finite area. As will be shown later, the extension over 
this sphere is possible. 



The case of a 



1/2 



In this case, we obtain 

g{t) = 1 + cxp (-V2<pt 



(37) 



The domain of time is — oo < i < oo. As in the case of 
Q 2 >l/2,i = — oo corresponds to an initial singularity at 
which g diverges. However, in contrast with the case of 
a 2 > 1/2, g is always larger than unity and finite except 
on t = — oo. This chart is incxtcndiblc. 



C. The case of a 2 < 1/2 
In this case, we obtain 



9( 



(i) = 1 + \/l - 2a 2 cosh (y/2ipt 



(38) 



The domain of time is also — oo < i < oo. As in the 
previous cases of a 2 > 1/2, t — > — oo corresponds to 
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an initial singularity at which g diverges. However, in 
contrast with the case of a 2 = 1/2, g diverges also in 
the limit of t — > oo. As will be shown later, the limit of 
t — > oo is not infinity, and thus this is a final singularity 
(big crunch). This chart is also inextendible. 



V. SOME PROPERTIES OF SOLUTIONS 

In this section, we shall analyze a few properties of this 
solution, including its cxtcndibility. For this purpose, we 
introduce a new time coordinate r defined by 



dr 



dt 



(39) 



The line element with this new time coordinate is then 
given by 



ds 2 = —(ipar) 2 dT 2 + dr 2 



(or) 



■dn 2 



(40) 



This new time coordinate is proportional to the proper 
time for a comoving observer along r =constant line. 
Thus, this is a geometrically and physically meaningful 
quantity. We discuss the cases of a 2 > 1/2, a 2 = 1/2, 
and a 2 < 1/2, separately. 



A. The case of a 2 > 1/2 

By integration of Eq. (|39p with an appropriate integra- 
tion constant, we have 



1 



2aip 



In 



+ yJ(V2a + l)/(V2a-l) 



1 - J(V2a - l)/(V2a + l)e^* 



We can easily see from the above equation that r vanishes 
for t — > — oo. This means that t — > — oo is not infinity. 
As mentioned, since the energy density p and Ricci scalar 
diverge for t — > — oo, this is the initial singularity By 
contrast, r becomes infinite in the limit of t — > As 
will be shown later, t = t\> corresponds to the future null 
infinity. 



B. The case of a" 



1/2 



In this case, we obtain 



C. The case of a 2 < 1/2 



In this case, we have 



2atp 



In 



1 + ^(1 + V2a)/(1 - V2a)e^« 



1 



(43) 



Also in this case, r vanishes in the limit of t — > — oo, and 
this "moment" corresponds to the initial singularity. In 
contrast to the case of a 2 = 1/2, we have, for t — ► oo, 



1 



2a(p 



In 



1 



1 - V2c 



(44) 



Thus, the limit of t — > oo does not corresponds to an 
infinity, but the final singularity. 
Using the above results, we obtain 

4a 2 e 2a V r 

q = — t= t= , (45) 

(V2a - l)e ia ^ T + 2e 2a ^ T - V2a - 1 

for arbitrary positive a. If the denominator in the right 
hand side of Eq. ([43]) vanishes, then g diverges. The 
denominator vanishes if 



1 



a 2aipr _ 



1 + s/2a 
1 - V2a 



(46) 



(47) 



is satisfied. The second root is meaningful only if a 2 < 
1/2, since r should be a real number. The first root 
is r = which corresponds to the big bang, whereas the 
second one is r = r c which corresponds to the big crunch. 
It might be a remarkable fact that, in the case of a 2 > 



(41) 1/2, the areal radius R behaves as 



i?-' 



V2c 



Aa 2 



■-r 2 e 2a *\ 



(48) 



for t ^> (atp)^ 1 . Thus it is seen that the comoving vol- 
ume exponentially expands as in the inflationary period; 
and that the acceleration of the cosmic volume expan- 
sion is realized, even though the strong energy condition 
holds. 

The world interval of the submanifold (r, r) in Eq. (|4"0|) 
takes the same form as that of the Rindler spacetime. 
Since the Rindler spacetime is extendible, we are led to 
infer that this solution might also be extendible. We 
discuss the maximal extension of this solution in the fol- 
lowing two sections. 



V2ip 



In 



,V2tpt 



(42) 



We can easily see that r vanishes in the limit of t — ► 
— oo, whereas r becomes infinite in the limit of t —> oo. 
Thus, t — > — oo is the initial singularity also in this case, 
whereas, as will be shown later, t — > oo corresponds to 
the future timclikc infinity. 



VI. ANALYTIC EXTENSION 

Following the prescription of the maximal extension 
for the Rindler spacetime, we introduce following new 
coordinates 

T = r sinh((/?ar) and X = r cosh((par) , (49) 
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then wc have 

ds 2 = -dT 2 + dX 2 + R 2 dn 2 . (50) 
From Eq. (|49[) . we have 

T 

X 2 -T 2 = r 2 and — = tanh(^ar). (51) 

X 

The first equation implies that r-constant curves are 
timelike hyperbolic curves, whereas r-constant curves are 
spacelike straight lines, in the (T,X) plane. The singu- 
larity r = corresponds to T = 0, whereas r = is null 
T = ±X, and thus T/X = ±1. This implies that r = 
also corresponds to r = ±co except at T = X = 0. 

The square of the areal radius R 2 is written as a func- 
tion of T and X in the form 

R 2 = T(V2aX-T). (52) 




t =const. 



r=0 (r=0) 



FIG. 1: The singularities and Hubble horizon of the analyti- 
cally extended solution are depicted in (T, X)-plane. 



The above equation explicitly shows the self-similarity of 
the present solution, since all the dimensionless variables 
are written as functions of a self-similar variable £ := 
X/T. Further, there is only one parameter a. Thus, the 
solutions obtained here compose a one-parameter family 
of self-similar solutions which is known as the Roberts 
solution (l2|. 

In order that R 2 is non-negative, both of T > and 
T < \J~2uX have to be satisfied, or both of T < and 
T > \/2aX have to be satisfied. The boundary of these 
two regions in (T, X) plane is the singularity T = (r = 
0). Thus we cannot regard these regions in (T, X) as two 
regions of one spacetime manifold: each region should 
be regarded as an independent spacetime manifold. Here 
the former, T > and T < V?aX, is of our interest. 

The energy density p is written in the form 



a 2 (X 2 



T 2 



2T 2 (V2aX - Tf 



(53) 







Thus the spacetime singularities arc located at T 
and T = \/2aX. 

It is useful for understanding the physical situation 
of this spacetime to examine the expansions of future- 
directed null, and also to consider the Misner-Sharp 
mass [l^ (see Appendix [Bj . In the coordinate system 
([50]) . the expansions of outgoing and ingoing null are 
given, respectively, by 



aX 
aX 



2 - a)T 
2 + a)T 



(54) 
(55) 



In accordance with Hayward's definition [40j, a region or 
a surface is said to be trapped if > 0, marginal if 

= 0, untrapped if ■d+d- < 0. A trapped region or 
a trapped surface of ??_ > is said to be past trapped, 
whereas a trapped region or a trapped surface of #_|_ < 
is said to be future trapped. Similarly, a region or surface 



of #_ = is said to be past marginal, whereas a region 
or surface of i? + = is said to be future marginal. 

From Eq. (|B2[) . we have the Misner-Sharp mass Mms 
of Eq. (jnnj) as 



M- 



1 MS 



(56) 



It should be noted that the Misner-Sharp mass is given 
by M MS = F/2, where F has been defined by Eq.Hr 



A. The case of a 2 > 1/2 

The domain covered by the original chart is < T < 
X. It is easily seen from Eqs. ([22]) and that T = X 
is regular (see also Appendix [Cj . Thus we can analyti- 
cally extend the region < T < X covered by the original 
chart (|3"2"j) over the region < T < \[2aX where the pos- 
itivity of R 2 is guaranteed. This extension have also been 
discussed in Oshiro et al [l4[ . R 2 vanishes on 
Here it should be noted that T = \/2aX is a spacetime 
singularity where the energy density p and scalar polyno- 
mials diverge there: see Eq. (f5"3"f and Appendix [Cl This 
singularity is a timelike naked singularity. 

It is seen from Eq. (|54|) that the expansion of out- 
going null is positive everywhere in the maximally 
extended chart, whereas the expansion of ingoing null 
i9_ is negative for aX/(y/2 + a) < T < s/2a, vanishes on 
T = aX/(V2 + a), and is positive < T < aX/(V2 + a). 
Thus, the region of < T < aX/(y/2 + a) is past 
trapped, T = aX/(y/2 + a) is a past marginal surface, 
and aXj (y/2 + a) <T < \/2aX is untrapped. The past 
marginal surface T = aX/ (v2 + a) corresponds to the 
Hubble horizon. The singularities and the Hubble hori- 
zon in (T, A)-plane are depicted in Fig.[TJ while the con- 
formal diagram is given in Fig. O Here it should be noted 
that the Hubble horizon is spacclikc; this is an important 



FIG. 2: The conformal diagram of analytically extended 
solution is depicted. This is similar to the upper half of 
the Minkowski spacetime. Note that the comoving lines, 
r =constant, enter into the future null infinity. This means 
that the elements of stiff fluid are accelerated to the speed of 
light asymptotically. 



difference between the present solution and inflationary 
solutions. 

We can easily see from Eq. (|5rj|) that Mms is positive 
for < T < X, Mms vanishes at T = X (r = in original 
coordinate chart), and Mms is negative for T > X. Fur- 
ther, Mms = ~ oo at R = 0, i.e., T — > y/2aX+o- The con- 
tribution of the scalar fields to Mms is positive and thus 
Mms is a monotonically increasing function with respect 
to X. Thus the negativity of Misner-Sharp mass in the 
region added by the extension comes from the negative 
infinite Misner-Sharp mass concentrated on the central 
timelike singularity. It seems difficult to get physical in- 
terpretation of this analytically extended solution, since 
the effect of the central naked singularity on its surround- 
ing region must not be negligible. 

The scalar field now takes the following form, 



V2 



T\/2a 2 - 1 



V2aX-T 



(57) 



At the singularities T = and T = s/2aX, the scalar 
field diverges. Here, it is interesting to see the norm of 
the gradient of the scalar field. We have 



g a %a^, b = ^(T 2 -X 2 ). 



FIG. 3: The singularities and Hubble horizon for a 2 = 1/2 are 
depicted in (T, X)-plane. There is a central null singularity 
which is not naked. 



V 



W7 



r=const. 

x + 

v z =const. 




Singularity'; 

T=o("=o) 



FIG. 4: The conformal diagram of analytically extended 
solution of a 2 = 1/2 is depicted. The comoving lines, 
r =constant, enter into the timelike infinity. 



B. The case of a 2 = 1/2 

In this case, from Eq. (|53|) . we find that the energy 
density p diverges at T = X, and thus T = X is a null 
spacetime singularity which is not naked, and its Misner- 
Sharp mass vanishes. As a result, the original chart is in- 
cxtcndiblc. Also in this case, T = A/3 is a past marginal 
surface which corresponds to the Hubble horizon. The 
singularities and the Hubble horizon in (T, A)-plane are 
depicted in Fig. 3. The conformal diagram is given in 
Fig. 4. By contrast to the case of a 2 > 1/2, comoving 
lines of constant r enter into a future timelike infinity. 
There is a future null and future timclike and spacelikc 
(58) infinities. 



We see from the above equation that </> ;a is timelike in 
the domain < T < X, null on T = X and spacelike in 
X < T < \/2aX. Therefore, this analytically extended 
solution is not equivalent to the system of a stiff perfect 
fluid. 



C. The case of a 2 < 1/2 

As mentioned earlier, in this case, there are two kinds 
of singularities in the universe. Both of them arc space- 
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o 



is? ,*' / ,\ a \N°v-' .' r= r=const. 



['_'- iBig bang singularity ^, 

r=o (t=0) 



FIG. 5: The singularities and Hubble horizon for < a < 
1/2 are depicted in (T, X)-plane. 



the extended region T > X can be Minkowski spacetime 
where the Misner-Sharp mass vanishes. In order to make 
clear whether such extension is possible, we consider the 
following Lorentz boost in the subspace (T,X), 



T=™ and X X + VT 



where 



v/2 

v = 1. 

a 

Then, we have the metric, for T — X < 0, as 
ds 2 = -df 2 + dX 2 + R 2 dn 2 , 

where 



(59) 



(60) 



(61) 



R 2 = X 2 - — (X - T) 



=const 




Singularity! 

FIG. 6: The conformal diagram of analytically extended solu- 
tion of a 2 < 1/2 is depicted. There are both the big bang and 
big crunch singularities although the topology of t =constant 
hypersurface is R 3 . 



like: one is the big bang and the other is the big crunch. 
It is remarkable that even though this is not a closed 
universe in the usual sense, the big crunch exists. There 
are the past and future marginal surfaces. Both of them 
correspond to the Hubble horizon, but the past one is 
in the expanding phase, whereas the future one is in the 
contracting phase. There is no null infinity. This type 
of solution is called the "universal" black hole by Carr 
and Gundlach 4l| . The same solution as this solution 
has also been obtained by Carr, Harada and Maeda[4^|. 
The singularities and past marginal surface are depicted 
in (T, A)-plane in Fig. 5. The conformal diagram is also 
given in Fig. 6. 



VII. C 1 EXTENSION FOR a 2 > 1/2: SPHERICAL 
VOID SOLUTION 



Here, we consider another extension for the case of 
a 2 > 1/2 without any serious singularities except for 
the big bang. For this purpose, it is very important to 
note that the Misner-Sharp mass vanishes on T — X = 0, 
i.e., (r, r) = (oo,0) in original chart. This implies that 



(62) 



By the Lorentz boost (|59p . the null hypersurface T— X = 
is mapped to T — X = 0, and thus we have R 2 = X 2 
at (r, r) = (oo,0). Further, we have, at (r, r) = (oo,0), 
or equivalcntly T — X = 0, 



drpR 2 = and d x R 



2X 



(63) 



Thus, if we put, for T + X < 0, the Minkowski spacetime 

ds 2 = -df 2 + dX 2 + X 2 dfl 2 , (64) 

we can easily see that this is a C 1 extension of the orig- 
inal solution. Note that the Ricci tensor of this solution 
does not vanish on T = X and hence the second order 
derivatives of the metric tensor are discontinuous. Also 
this extension have been discussed in Oshiro et al [lij]. 

The region newly added by this extension is regarded 
as an expanding void. The big bang singularity and cos- 
mological horizon of this C 1 extended spacetime are de- 
picted in (T, A)-planc in Fig[7J The conformal diagram 
of this C 1 extended spacetime is depicted in Fig[Hl 

This universe model has a few very peculiar properties. 
Since the gradient of the scalar field is null on T = X, 
this solution cannot be regarded as the universe filled 
with a stiff perfect fluid only. However, we can regard 
the matter field in T < X as the stiff perfect fluid. The 
world lines of fluid elements are the curves with constant 
r. These are expressed in (T, A)-coordinate system by 
hyperbolic curves 



X 2 -T 2 = r 2 . 



(65) 



We see from the above equation that the world lines 
of fluid elements become asymptotically null. This fact 
might be rather surprising, since the fluid elements are 
accelerated outward even though there is a positive grav- 
itational mass (Misner-Sharp mass) inside the mass shell 
on which the fluid elements stay. 
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r=const. 
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FIG. 7: In the case of C 1 -extended solution, the big bang 
singularity, cosmological horizon and light ray are depicted in 
(T, X)-plane. 



Newly Added 
Minkowski Region 



r=const. 



Original Chart 




r=const. 



Singularity t=^>X(t=0) 



FIG. 8: The conformal diagram of C 1 maximally extended 
solution is depicted. This solution describes an expanding 
spherical void in the universe filled with a massless scalar- 
field. 



The outward acceleration causes the temporal growth 
of the redshift of the light ray emitted from a fluid ele- 
ment in T < X to the observer at the symmetry center 
X = 0. (see Fig J51) The ingoing null geodesies are given 
by 



T = ujqX + To and X 



-uj X, 



(66) 



where A is the affine parameter, ujq and To are positive 
constants which correspond to the angular frequency de- 
tected by the observer at (T, X) = (T , 0). The big bang 
singularity is given by T = —vX in this coordinate sys- 
tem, and thus the affine parameter takes the values in 
—To/y/2ujo(l — v) < A < 0. The tangent vector of this 
null geodesic is 



k» = (w , -w 0) 0,0), 



(67) 



The 4-velocity of the fluid element labeled by r is given 
by 



" 2 + r 2 ,T,0,0 



(68) 



The light ray detected by the observer at (T, X) = (T , 0) 
should be emitted from the fluid element of r at 



(T,X) = (f e ,X e ) := 



T 2 



T 2 



2Tn 



2T n 



Thus the redshift z of this light ray is given by 



~k a " u a\(T,X) = (T <! ,X <! 
U) 



1 



(69) 



(70) 



The important and peculiar feature is that the redshift 
z is monotonically decreasing function of r. For a fluid 
element of To < r, the redshift z is negative, i.e., the light 
suffers the blucshift. The fluid element of To < r emits 
the light ray at T = T c < 0, and at this moment, the 
radial component of the 4-velocity is negative. Thus the 
reason of the blucshift for To < r is recognized as a result 
of the ingoing motion of the fluid element. By contrast, 
if To > r, i.e., T c > 0, the redshift z is positive. If To = r, 
i.e., T e = 0, the redshift vanishes. The temporal variation 
rate of z is 



# = i>o. 

dT Q r 



(71) 



The redshift temporally increases. This is also a very dif- 
ferent property from the dust-filled inhomogeneous uni- 
verse mo del [H, HI]. 

Since the cosmological horizon T = (1 — l/a 2 )X, or 
equivalently, T = aX/{-j2 + a), is spacelike, the fluid el- 
ements necessarily enter into the inside of the cosmolog- 
ical horizon from its outside. Thus even though the fluid 
elements are accelerated outward, the inflation does not 
occur in usual sense. (We note that in the inflationary 
period, comoving world lines go outside the cosmological 
horizon from its inside, and such behaviors of comoving 
world lines are essential to resolve the horizon problem.) 

Finally, we investigate the distance-redshift relation. 
Observationally useful distance is the luminosity distance 
c?l[H which is given by 

rf L = (1 + z) 2 R\(f,X) = (f ,X ) 

= ^V(l + z) 4 + 2(l + z) 2 -2a 2 + l. (72) 
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Since c?l is monotonically increasing function of z, it is 
monotonically decreasing function with respect to r. The 
redshift of the light ray from the big bang z = z\> is given 
by 

z h = yfV2a-l-l. (73) 

We can easily see that the luminosity distance of the 
big bang vanishes. Since the light ray emitted from a 
comoving source at the big bang suffers finite blucshift. 
the flux from the big bang Tb = Lb/ind^ diverges at 
the observer, if the luminosity of the big bang singu- 
larity Lb is finite. This feature implies that comoving 
sources of radiation behave as white holes at this big 
bang sing ularityjH, IH|. In the limit r — > 0+, z — > oo 
and thus c?l of the boundary of void r = is infinite. 
It is hard to observe the vicinity of the boundary of the 
void, while it is easy to observe the vicinity of the big 
bang, in this universe. 

VIII. SUMMARY AND DISCUSSION 

We have rederived and investigated here the solution 
obtained by Roberts, which is a spherically symmetric 
but inhomogencous universe filled with a massless scalar 
field minimally coupled to gravity, from a cosmologi- 
cal perspective. The solutions obtained compose a one- 
parameter family, which is divided into three distinct 
classes. The first class consists of solutions with only one 
spacelike singularity in the comoving chart. The second 
class consists of solutions with two singularities which 
are null and spacelike, respectively. The third class con- 
sists of solutions with two spacelike singularities which 
correspond to the big bang and big crunch, respectively. 

In the case of the first class, the comoving chart does 
not cover the whole spacetime. Hence, we constructed 
two maximally extended solutions from the solution in 
this class. The analytic extension leads to the solution 
which contains a timclikc singularity at the symmetry 
center and thus seems to be unphysical. By another ex- 
tension, we obtained a solution which has no singularity 
other than the big bang but one that contains a spheri- 
cal void. The later one has very peculiar but interesting 
properties. If the gradient of the massless scalar field is 
timelike, this is equivalent to the stiff perfect fluid case. 
Then in the region where the scalar field has a timelike 
norm, we can naturally define fluid elements and consider 
their motions. In this solution, the fluid elements move 
outwards, and further, their outward speeds are accel- 
erated, even though the strong energy condition holds. 
This is a feature which is significantly different from the 
homogeneous and isotropic universe models filled with 
the matter which satisfies the strong energy condition, 
and also from spherically symmetric but inhomogeneous 
universe filled with the dust matter. 

The physical reason of the outward acceleration of the 
universe that we have deduced here appears to lie essen- 
tially in the inherent physical nature itself of the massless 



scalar field. To try to understand this better, we note 
that we can construct a similar and somewhat parallel 
situation in a Minkowski background. Using the same 
coordinate system as Eq. , let us consider a following 
solution, 

0(T,X) = i(X-T) 2 H (X-T) (74) 

where A is constant and 6*h is the Heaviside's step func- 
tion. We focus on the domain of T > 0, since <j> is every- 
where finite for T > 0, but is infinite at X — for T < 0. 
Then we have 

g ab 4>-A-,b = - T) 3 (3A + f)9 H (X - f). (75) 

We can easily see from the above equation that <p. a is 
timelike for T < X, whereas it is null at T = X. Further, 
we see that in the limit of T — > oo with U := T — X fixed, 
the norm of </> ;a becomes 

9 ab 4>; a <t>; b -> jim^ A \ U 3 (4f + 3U) = 0. (76) 

The above results shows that (f>- a becomes null asymp- 
totically. This behavior is very similar to the present 
solution. This seems to indicate that the outward accel- 
eration comes from the inherent nature of the massless 
scalar field. As it turns out, in the present solution, the 
gravity produced by the scalar field itself is too small to 
decelerate the outward motion of the scalar field, even 
though the total (Misncr-Sharp) mass for the universe is 
infinite. 
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APPENDIX A: ON THE FUNCTION f(r) 

For notational simplicity, we define the following quan- 
tities 

X-=r 2 f{r) and it := g(t) +p(r), (Al) 

where / and p are defined by Eqs. (|19|) and respec- 
tively. Then, Eq. (f!?5|) is written in the form 



H + l = 
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By differentiating the above equation with respect to r, 
we have 



CqU 

where 



c 5 u 



c 3 u 



Cl 
C2 

C3 

f'4 

C5 

C6 



-4/i'V, 
2 (h'h" - 



= 2 



hp' 

X 



c 2 u + ci = 0, 



3xW) 

-x P 



p- 



(A3) 

(A4) 
(A5) 

(A6) 

(A7) 

(A8) 
(A9) 



Here note that all of the coefficients c„ do not depend 
on t. This implies that if there is a non-trivial real root 
u = U of Eq. (|A3[) . U must depend on only r. If so, we 
have 



g(t) = U(r)-p{r). 



(A10) 



The above equation holds only if both of g and U — p 
are constants. Thus, there are two possibilities; one is 
g(t) ^constant, and the other is that all of the coefficients 
c n (n = 1,..,6) identically vanish so that Eq. (|A3[) is 
trivial. In the latter case, g may depend on t. 

In the former case, without loss of generality, we can 
put u = p. Then, by definition of p, we find 



ti 



(All) 



Further, substituting g = into Eq. (|27[) . we have H = 0. 
By using this fact and Eq. (|A11|) . x can be eliminated 
from Eq. (|A3|) . and we have 



From C4 = 0, we have x' ^constant, and thus 

X = ar + /?, (A15) 
where a and (3 are integration constants. Eqs. (|A14[) 
and (|A15[) lead to c\ = c 2 = C3 = 0. Since we assume 
that r = is a non-singular point at least initially, R 
vanishes at r = 0. Thus, we can see from Eq. (jTSJ) that 
(3 should vanish. Further, it is seen from Eq. (| 14[) that F 
should vanish at the regular origin r = 0. Hence, from 
Eq. (f2"Tj) . h should vanish at r = 0, and thus, we have 
/i = 0. 



APPENDIX B: MISNER-SHARP MASS 

In general, the metric of the spherically symmetric 
spacetime is given by 

ds 2 = A 2 {t,r)dt 2 - B 2 {t,r)dr 2 - R 2 {t,r)dn 2 . (Bl) 

In this coordinate system, the Misner-Sharp mass is given 
by[li 



7? / R 
M MS = - I 1 + — 0+0- 



where 



(B2) 



(B3) 



Here note that 1?+ is the expansion of the outgoing radial 
null, while i?_ is the expansion of the ingoing null. 



APPENDIX C: RICCI TENSOR 



P 



h" , 



P 



1 /4 
P P 



1 



hp , 



El 

2 



h>' v j h' 

(A12) 

The above differential equation gives a relation between 
h(r) andp(r). The above equation, however, corresponds 
to the static case, in which we are not interested here in 
the present work. The static solutions will be discussed 
elsewhere. 

As for the dynamical case, the conditions c„ = give 
differential equations for x and h. From cq = 0, we have 



h , 

— X = °- 

X 



Thus, we have h - 
hh' = 0, and thus 



(A13) 

or x' — 0. From C5 = 0, we have 



constant. 



(A14) 



Non-vanishing components of the Ricci tensors with 
respect to the analytically extended chart are 



d d 

"[dT'dX 



a 2 TX 



Ric 



R 4 

d d \ a 2 T 2 



dX'dXJ i? 4 
The Ricci scalar is then given by 



trffic=— (T 2 -X 2 ). 



(Cl) 
(C2) 
(C3) 

(C4) 



Since R vanishes on T = \/2aX , the Ricci scalar trffic 
diverges there. 
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